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ABSTRACT 

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of 
black holes and host galaxies predicts quasar luminosity functions at redshifts 1 -5 < z < 6 and local properties in 
nice agreement with observations. In our model the quasar activity marks and originates the transition between 
an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, 
and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while 
the latter accounts for the high number density of massive galaxies at substantial redshifts z > L5, the popula- 
tion of Extremely Red Objects, and the properties of local ellipticals. 

Subject headings: galaxies: formation - galaxies: evolution - quasars: general 



1. INTRODUCTION 

Observations in the optical and in the X-ray bands have shown since a couple of decades that very bright quasars (QSOs) occur 
very soon in the history of the Universe, and that their average luminosity significantly declines from redshift z ^ 2 toward the 
present (e.g., Mathez 1978; Giacconi 1985; Schmidt & Green 1986). Modern observations have confirmed and detailed such a 
picture (e.g.. Fan et al. 2004, 2006; Richards et al. 2005, 2006; Cristiani et al. 2004; Barger et al. 2005; Tozzi et al. 2001; Brandt 
& Hasinger 2005), measuring QSO luminosity functions (LFs) up to z ^ 6. 

Supermassive black holes (BHs) were found to be ubiquitous in the centers of spheroidal galaxies (Kormendy & Richstone 
1995). Moreover, a narrow relationship between the central supermassive BH mass and the stellar mass/luminosity of the 
spheroidal/bulge component has been firmly established by observations (e.g., Magorrian et al. 1998; Marconi & Hunt 2003). 
The BH mass also strictly correlates with the velocity dispersion of the spheroidal component (e.g., Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Tremaine et al. 2002). Contrariwise, it has been shown that the BH mass is on average at least 10^ times 
smaller in late-type irregular and/or spiral galaxies than in ellipticals with the same luminosity (see Salucci et al. 2000). All 
together these observations witness the strict connection between the mass of the BH and that of the old stars formed at z > 1 
within the host galaxy. Recently, the studies of QSO host galaxies have been extended to high redshift (e.g., Dunlop et al. 2003; 
Floyd et al. 2004), with the general result that they are typically early-types. 

The stellar populations of ellipticals are old and have similar ages/formation histories, although there is some trend for more 
massive systems to be on average older and more metal-rich (Sandage & Visvanathan 1978; Bernardi et al. 1998; Trager et 
al. 2000; Terlevich & Forbes 2002; Thomas et al. 2005). A color-magnitude relation is well established, indicating that more 
luminous spheroids are redder (Bower et al. 1992). The widely accepted interpretation is that brighter objects are richer in metals 
and that the spread of their star formation epochs is small enough to avoid smearing of their colors. In fact, the slope of the 
relation seems not to change with redshift (Ellis et al. 1997; Kodama et al. 1998), supporting this view. Note, however, that the 
trend of having on average older stellar populations in more massive ellipticals may also constitute an important ingredient (e.g., 
J0rgensen et al. 1999; Caldwell, Rose & Concannon 2003; Nelan et al. 2005). 

The star formation history of spheroidal galaxies is mirrored in the Fundamental Plane (e.g., Jedrzejewski et al. 1987; Dressier 
et al. 1987) and in its evolution with redshift. Ellipticals adhere to this plane with a surprisingly low orthogonal scatter (around 
15%), as expected for a homogeneous family. Furthermore, recent studies (e.g., Treu et al. 2002; van der Wei et al. 2004; Holden 
et al. 2005) suggest that ellipticals, both in the field and in the clusters, lay on the fundamental plane up to z ~ 1, consistent with 
massive spheroids being old and quiescent. Progenitors of local massive early-type galaxies have been identified through K-band 
and Spitzer surveys at substantial redshifts, z ^ 1 • Direct evidence that massive galaxies with stellar content M^, > 10" Mq were 
in place at z > 2 is provided by recent K-band surveys (e.g., Cimatti et al. 2002; Kashikawa et al. 2003; Fontana et al. 2004; 
Bundy, Ellis & Conselice 2005). The space density of Extremely Red Objects (EROs) at z > 3 is only a factor about 5-10 
less than that at z ^ 1 (e.g., Tecza et al. 2004). This implies a phase of extremely high star formation rate (SFR), from several 
hundreds to thousands solar masses per year, taking place at high z. Such a phase is also witnessed by the sub-mm galaxy counts 
(e.g.. Chapman et al. 2003; 2005). 

On the theoretical side, the hierarchical clustering paradigm led to the development of various semi-analytic models for galaxy 
formation. These generally share the basic assumption that gravity is the main driver in shaping the structure and morphology 
of galaxies. In this scenario, the gas cools and form stars following the collapse of dark matter (DM) halos. In the standard 
cosmology, small DM objects form first and merge together to make larger ones. This scenario then implies that large ellipticals 
form late, by mergers of disk/bulge systems made primarily of stars, recently indicated as 'dry mergers' (e.g., Naab et al. 2006). 
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But the early appearance of QSOs with high luminosity and presumably huge BH mass is at variance with expectations from 
this framework (e.g., Bromley et al. 2004). This apparent contradiction may be solved by invoking a higher efficiency in forming 
massive BHs in smaller galaxy halos at higher redshift (e.g., Haehnelt & Rees 1993; Haiman & Loeb 1998; Wyithe & Loeb 2003; 
Mahmood et al. 2005). 

Granato et al. (2001, 2004), instead, explored the possibility of reconciling the observed 'downsizing' of QSOs and spheroidal 
galaxies with the bottom-up DM hierarchy (see also Somerville et al. 2004; Baugh et al. 2005; Croton et al. 2006; Hopkins et 
al. 2006). This has been done by developing a simple model that incorporates the main physical aspects of the DM and baryons 
residing within galactic halos. This new approach emphasized the role of energy feedback processes both from supernovae 
(SNae), capable of unbinding the gas in low-mass systems, and from the QSO phase of the supermassive BHs, capable of 
ejecting gas from the largest objects (see Di Matteo, Springel & Hemquist 2005; Lapi, Cavaliere & Menci 2005). These feedback 
processes can actually reverse the formation sequence of visible galaxies with respect to that of DM halos. Hence large galaxies 
end their star formation, and have their BHs shining as QSOs early on. On the contrary, the star formation and the QSO phase 
are more prolonged in smaller halos (hence the name of 'Anti-hierarchical Baryon Collapse', or ABC, scenario). 

While in previous works we focused on reproducing the properties of spheroidal galaxies, in this paper we explore the con- 
straints on model parameters set by the redshift-dependent LFs of QSOs. Our plan is as follows: in § 2 we briefly recall the 
main features of the model by Granato et al. (2004); in § 3 we describe the procedure adopted to compute the QSO LFs, the 
supermassive BH mass function, and other galactic observables; § 4 is devoted to the illustration of our results; in § 5 we make a 
critical comparison of our findings with other models in the literature; finally, in § 6 we summarize and discuss our conclusions. 

Throughout this work we adopt the cosmology indicated by the WMAP data (Bennett et al. 2003; Spergel et al. 2006), i.e., a 
flat Universe with matter density ~ 0.27, baryon density fib « 0.044 and Hubble constant //q ~ 71 km s"' Mpc~'. 

2. THE MODEL FROM THE GROUND UP 

This paper is based on the semi-analytic model developed by Granato et al. (2004), which follows the evolution of baryons 
within proto-spheroids through simple and physically grounded recipes. 

We defer the interested reader to that paper for a full account of the physical justification and a detailed description of the 
model. Here we provide a short summary of its main features, focusing on the aspects relevant to our discussion on QSO LFs 
and BH demographics. 



2.1. The dark matter sector 

As for the treatment of the DM in galaxies, the model follows the standard hierarchical clustering framework, taking also into 
account the results by Wechsler et al. (2002), and Zhao et al. (2003a; 2003b). Their simulations have shown that the growth of a 
halo occurs in two different phases: a first regime of fast accretion in which the potential well is built up by the sudden mergers 
of many clumps with comparable masses; and a second regime of slow accretion in which mass is added in the outskirts of the 
halo, without affecting the central region where the galactic structure resides. 

This means that the halos harboring a massive elliptical galaxy once created, even at high redshift, are rarely destroyed. At 
low redshifts they are incorporated within groups and clusters of galaxies. Support to this view comes from studies of the mass 
structure of elliptical galaxies, which are found not to show strong signs of evolution since redshift z « 1 (Koopmans et al. 2006). 
Note that, as pointed out in § 1, the BH mass is strictly correlated with properties (mass and velocity dispersion) of the old stars 
in massive early-type galaxies, formed at least 8 Gyr ago (see Thomas et al. 2005) and, as a consequence, in massive galaxy 
halos virialized at z > 1 .5. 

We deal just with the latter redshift range (z > 1 .5), where a good approximation of the halo formation rates is provided by the 
positive term in the cosmic time derivative of the cosmological mass function (e.g., Haehnelt & Rees 1993; Sasaki 1994). For 
DM halos with mass Myk at time fyir, these formation rates are given by 



d^A^s 
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where NsriMyi,-,!) is the Sheth & Tormen (1999, 2002) version of the PS mass function (Press & Schecter 1974). In the above 
equation, a = 0.707 and p = 0.3 are constants obtained from comparison of the mass function with the outcome of A^-body 
simulations; c7(Mvir) is the mass variance of the primordial perturbation field, computed from the Bardeen et al. (1986) power 
spectrum with correction for baryons (Sugiyama 1995), and normalized to trg w 0.84 on a scale of Sh~^ Mpc; ddtm) is the 
critical threshold for collapse, extrapolated from the linear perturbation theory. Note that our adoption of the Sheth & Tormen 
mass function to construct the rates is mandatory, since it is well-known that the canonical Press & Schechter theory strongly 
underpredicts the number of massive halos, particularly at the high redshifts relevant for the computation of the QSO LFs. 

We confine our analysis to galaxy halo masses M™ ^ 2 x 10^ 'M©, since we are interested in following the evolutionary 
history of bright QSOs and of their host galaxies. At the other end, weak lensing observations (e.g., Kochanek & White 2001; 
Kleinheinrich et al. 2004) and kinematical measurements (e.g., Kronawitter et al. 2000; Gerhard et al. 2001) suggest an upper 
limit in galaxy halo mass M™" « 2 X IO'^Mq. The same limit is also indicated by the theoretical analysis of Cirasuolo et 
al. (2005) on the velocity dispersion function of early-type galaxies. These limits in mass and redshift ensure that the positive 
cosmic time derivative of the halo mass function is a good approximation to the formation rate of DM halos, as the negative term 
is negligible. 

Note, however, that since we are interested in DM halos associated to a single galaxy, we need to correct the Sheth & Tormen 
halo mass function by: (i) accounting for the possibility that a DM halo contains various sub-halos each hosting a galaxy; (ii) 
removing halos corresponding to galaxy systems rather than to individual galaxies. We deal with (ii) by simply cutting-off the 
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Fig. 1 . — Left panel: dependence on galactic age of the stellar mass, infalling mass, cold gas mass, and star formation rate within halos of various masses, 
virialized at redshift 4. Right panel: same for the BH mass, reservoir mass, and BH accretion rates. 



halo mass function at M™" w 2 x IO'-'Mo, beyond which the probability of having multiple galaxies within a halo quickly 
becomes very high (e.g., Magliocchetti & Porciani 2003). As for (i), we add the sub-halo mass function, following the procedure 
described Vale & Ostriker (2004; 2006) and Shankar et al. (2006), and using the fit to the sub-halo mass function at various 
redshifts provided by van den Bosch, Tormen & Giocoli (2005). However we have checked that for the masses and redshifts 
relevant here (z ^ 1 .5 - 6 and Mvk ^2x10" Mq - 2 x lO'^M©), the total mass function (sub-halo + halo) differs from the halo 
mass function by less than 5%. Thus, using simply the halo mass function to derive the DM halo formation rates is a good 
approximation. 

2.2. The baryonic sector 

The physics governing the evolution of the baryons is much more complex than for the DM. The main features of our model 
can be summarized as follows (see Granato et al. 2004 and Cirasuolo et al. 2005 for additional details). 

During or soon after the formation of the host DM halo, the baryons falling into the newly created potential well are shock- 
heated to the virial temperature. The hot gas is (moderately) clumpy and cools quickly especially in the denser central regions, 
triggering a strong burst of star formation. The radiation drag due to starlight acts on the gas clouds, reducing their angular 
momentum. As a consequence, a fraction of the cool gas can fall into a reservoir around the central supermassive BH, and 
eventually accretes onto it by viscous dissipation, powering the nuclear activity. The energy fed back to the gas by SN explosions 
and BH activity regulates the ongoing star formation and the BH growth. Eventually, most of the gas is unbound from the DM 
potential well, so that the star formation and the BH activity come to an end on a timescale shorter for the more massive galaxies. 

In Appendix A and Table Al we present the basic equations and parameters controlling the evolution of the baryonic component 
in our model, once the halo mass and the virialization redshift are given. These equations can be numerically integrated to yield, 
among others, the SFR and the accretion rate onto the central BH as function of cosmic time. In Fig.^we plot the basic outputs 
of the model. 

Initially, the cooling is rapid and the star formation is very high; thus the radiation drag is efficient in accumulating mass into 
the reservoir. The BH starts growing from an initial seed with mass 10^ M0 already in place at the galactic center. Since there 
is plenty of material in this phase, the accretion is Eddington (or moderately super-Eddington) limited (e.g.. Small & Blandford 
1992; Blandford 2004). This regime goes on until the energy feedback from the BH is strong enough to unbind the gas from 
the potential well, a condition occurring around the peak of the accretion curve. Subsequently, the SFR drops substantially, the 
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radiation drag becomes inefficient, the storage of matter in the reservoir and the accretion onto the BH decrease by a large factor. 
The drop is very pronounced for massive halos Mvii > 10^^ Mq, while for smaller masses a smoother declining phase can continue 
for several Gyrs, and the BH and stellar masses can further increase by a substantial factor. 

Before the peak, radiation is highly obscured by the surrounding dust. In fact, these proto-galaxies are extremely faint in the 
UV-optical rest frame and are more easily selected at (sub-)mm wavelengths. The nuclear emission is heavily obscured too; 
however, since the absorption significantly decreases with increasing X-ray energy of photons, it is easier to detect it in hard 
X-ray bands. On the other hand, in the proximity of the accretion peak, i.e., when the central supermassive BH is massive and 
powerful enough to remove most of the gas and dust from the surroundings, the active nucleus shines as an optical QSO. 

The list of the observations well fitted by the model is presented in Table A2, along with a brief description of the 
main underlying assumptions. In order to translate the SFR and the mass in stars into more directly observable quan- 
tities such as broad-band luminosities, including dust extinction, we exploit the GRASIL code (Silva et al. 1998, see 
jhttp : //web .pd. astro . it/granatol. As for the QSOs, the model yields bolometric luminosities that are translated 
into luminosities in a given band by means of appropriate bolometric corrections. The epoch-dependent luminosity functions are 
derived combining the results obtained for halos of fixed mass and virialization redshift with the halo formation rates presented 
in §2.1. 

3. THE MODEL AT WORK 

As a first step we want to properly translate the DM formation rates given by Eq. Q into BH formation rates. This calls for a 
relation between the halo mass Mvir and the BH mass M,. In Fig.|2](left panel) we plot the one obtained by solving the equations 
presented in Appendix A, that can be well approximated by 

^ (M..,/2.2xlO"Mof'^^ 

This relationship is very close to that found, at z = 0, by Shankar et al. (2006) by comparing the local BH mass function with the 
galaxy halo mass function. However, this is only an average relationship and we do expect a significant scatter around it, since 
the values of the parameters in the equations of Appendix A may naturally vary from halo to halo. In the following we assume 
that the above relation logM, -logMvir holds on average with a gaussian dispersion AlogM,. Therefore we convert the halo 
formation rates into BH formation rates through the convolution 

d2 -(logM;-logM.)-/2(AlogM.)^ 

. (3) 
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It will be shown below that the dispersion is an important ingredient when comparing model predictions to observations. 

Note that despite the redshift dependence of the M,/Mvii ratio (see Eq. (21), our model predicts a relatively weak evolution of 
the BH-bulge mass ratio. This is because, as mentioned in § 2.2, the growth of the BH mass is controlled by the SFR through the 
radiation drag and the SN feedback; in turn, the feedback from the QSO can eventually sweep out the gas halting both the star 
formation and the accretion on the BH, especially in massive galaxies. In the latter structures the stellar and BH masses grow 
(and stop growing) in parallel, while in smaller objects where only SN feedback is effective the BH-bulge mass ratio turns out to 
depend slightly on z (see Fig. 5 of Cirasuolo et al. 2005). 



3.1. QSO luminosity functions and BH mass function 
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The QSO LFs can now be computed. Up to its peak, the BH bolometric light curve can be well approximated by the simple 
exponential form 



AM 

m = — e('-'--^'P-)/-- Oh (fv„ + Afpeak - Af^.s < t < + Afpeak) • (4) 

fEdd 

Here fEdd « 4 x 10*^ yr is the Eddington timescale, and Tef « ?7fEdd/(l - 77) A is the e-folding time in terms of the BH mass-energy 
conversion efficiency and of the Eddington ratio A. The Heaviside function^ 6^ specifies that the QSO shines unobscured only 
during the time interval Afyis before the peak of its light curve. This means that, as discussed in § 2.2, at the time rvir+ Af peak- Af vis 
the BH is massive enough for its energy feedback to remove most of surrounding gas and dust, letting the active nucleus to be 
visible. In this work we have used an Eddington ratio A slightly rising toward high z, and specifically A = 4 for z > 6, A = 3 for 
5 < z < 6, A = 1 .7 for 3 < z < 5, A = 1 for 2 < z < 3, and A = 0.8 for 1 .5 < z < 2. The empirical formula A(z) « -1 . 15 + 0.75 (1 +z) 
works as well in the redshift range 1 .5< z < 6 (see also Appendix A and Table Al). This assumption is discussed in § 6. 

There are three timescales in Eq. Q: the virial time fvii depends on cosmology; the peak time Afpeak is obtained by solving 
the system of equations reported in Appendix A; the visibility time Afvis, dependent on dust and gas absorption, is taken as a 
parameter, since its computation is challenging for semi-analytical models. 

The QSO LF at a time t and luminosity L is computed by summing up the contributions of aU the sources which virialize at 
epochs fvii ^ t and shine at the time t with luminosity L. One has 

J,-Ar,^^ J dfvir dM, V fEdd / 

where 6d indicates the Dirac delta function. The time delay Afpeak provided by our code ranges from 0.2 Gyr at redshift z > 5 
where A « 4 to 1 Gyr at redshifts z < 2, where A < 1 (see Fig.|2l right panel). 

In our model, the declining phase of the QSO light curve can be neglected for halo masses above few x IQ^^Mq, because for 
them the QSO feedback is so powerful to remove suddenly most of the accreting gas from the surroundings. As a consequence, 
after the peak the accretion rate drops abruptly to values much lower (by factors of 10"^ at least; see Fig. [3 than the Eddington 
rate, and the radiation process becomes much less efficient (for a review see Blandford 2004). We have checked that such 
radiatively inefficient stages provide a negligible contribution to the QSO LFs, in accord with the findings by Yu & Lu (2004) 
based on the comparison between the observed QSO LFs and the local supermassive BH mass function. 

On the contrary, for halo masses smaller than 1O'^M0 and especially at redshifts z < 2, the BH accretion rate decays more 
gently after the peak, and up to few Gyrs are required before it lowers to levels of ~ 10~^ the Eddington rate, becoming radiatively 
inefficient. This declining phase is consistent with observations of moderate QSO activity in relatively typical early-type galaxies, 
with old stellar populations (e.g., Kauffmann et al. 2003). This regime, which adds to the visibility time defined above (that refers 
to the rising phase of the bolometric light curve), is included in our computation of the QSO LFs. We find that it contributes 
significantly to the faint end of the X-ray QSO LFs (see also Hopkins et al. 2006b) and to the X-ray number counts (see Granato 
et al. 2006 for details). On the other hand, the optical luminosities are generally below those corresponding to the standard 
definition of optical QSOs (Mg < -23); therefore the declining phase has a lower effect on optical LFs. 

The last step of our computation is the conversion of the bolometric QSO LFs computed above to the optical and hard X-ray 
bands through the appropriate bolometric corrections (e.g., Elvis et al. 1994; Ueda et al. 2003; Vignali, Brandt & Schneider 
2003; McLure & Dunlop 2004; Barger et al. 2005; Richards et al. 2005, 2006). For the optical band we use L/Lb = fs with 
fs ~ 10. For the hard X-ray band (2-10 keV) we take into account the observed luminosity dependence of L/Lx = fx, given by 
fx ~ k]l^^^^^^\L/m^^txg s-i)/3''/(ft+i), with kx ~ 17 and Px ~ 0.43 (see also Shankaret al. 2004). 

The epoch-dependent BH mass function is given by 

mogM„t)= / dfv„- . . (6) 

7o dfvii-dlogM, 

3.2. Related galactic properties 

Our model aims at predicting, at the same time, both the QSO LFs at various redshifts, and the properties of host galaxies. 
Using the DM halo formation rates of Eq. Q we can compute the galaxy LF in a given band as 



d^A^cT 
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dfvir dMv 
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where Li is the luminosity in the i-th band. Similarly, the velocity distribution function is given by 

on A Th. '^'^st dMv 
iZ(logcr,f)= / df 



dfvir dMv 



dlogcr 



(7) 



(8) 



To derive the Jacobians we use: (i) the GRASIL code, which evaluates the luminosity of a galaxy as a function of time once 
the dust content and the SFR (obtained by solving the system of equations in Appendix A) is specified; (ii) the scaling law for 

^ The Heaviside function 9h is defined by 

1, if a: is true; 



0, otherwise . 
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Fig. 3. — Left panel: predicted QSO LF in the optical B-band at z = 6 (solid line) and z = 8 (dotted line) for an Eddington ratio A = 4, and a visibility time 
Afyis Ki5 X 10^ yr. Data points are from Fan et al. (2004). Right panel: QSO LF in the optical B-band at z = 5 for A = 3 and A?vis ~ 5 X 10^ yr. Data points are 
from Kennefick et al. (1995, open squares); Fan et al. (2001, diamonds and asterisks); Wolf et al. (2003, triangles); Cristiani et al. (2004, circles); Richards et 
al. (2006, filled squares). 



the virial velocity Vvii- = (GMyh/Z^vir)'''^ cx m\[^ (1 +Zvir)''^^; (iii) the ratio between virial velocity and stellar velocity dispersion, 
a = 0.55 Vvir, found by Cirasuolo et al. (2005). The combination of these three ingredients yields a magnitude-velocity dispersion 
(Faber-Jackson) relation in good agreement with that (M, * = -7.5 log (T-5) observationally derived by Sheth et al. (2003) through 
the bisector method (see also Cirasuolo et al. 2005). 

Exploiting the $(Li,f), it is straightforward to compute the galaxy counts and the related redshift distributions. These estimates 
have been worked out in particular for sub-mm and K-band selected galaxies, whose counts and redshift distributions proved to 
be particularly challenging for competing semi-analytic models. We defer the interested reader to the papers by Granato et al. 
(2004) and Silva et al. (2005) for the description of the computational strategy, based on the GRASIL code, and of the results. 

4. RESULTS 
4. 1 . QSO Luminosity Functions 

An important feature of the model is the prediction of a pre-QSO phase lasting Afpeak ^ 15-20e-folding times (for the most 
massive objects), during which the mass of the central BH is increasing exponentially. Meanwhile the host galaxy is forming 
stars at extremely high rates (cf. Appendix A and Fig.^ see also Sanders et al. 1988; Croom et al. 2006; Hopkins et al. 2006) 
and is therefore very bright at far-IR and (sub-)mm wavelengths. 

In Figs. l3l6l we show that our model provides a very good fit to the optical and hard X-ray QSO LFs over the full redshift range 
over which it applies (z > 1.5). Due to the delay Afpeak between the halo virialization epoch and the peak of the BH accretion, 
the QSOs shining at the cosmic time tniz) are associated to halos virializing at an earlier time fvir ~ fH(z)- Afpeak < tniz)- 

This is an important point since, at high redshifts and for massive halos, Afpeak corresponds to a substantial fraction of the 
cosmic time (see Fig. |2] right panel), and the DM halo formation rates at fvii are significantly lower than at tniz)- This is 
illustrated by Fig. |3l (left panel), where we also show the predicted optical LF of QSOs at z = 8. The precipitous drop of the 
number density of bright QSOs, compared to z = 6, is due to the fact that, for tniz = 8), Afpeak ~ 0.2 Gyr corresponds to Az « 2, 
so that the abundance of bright QSOs reflects that of virialized massive halos at z ~ 10, which is far lower than at z ~ 8. 

The flattening of the LFs at the low end is due to the flatter slope of the formation rate of less massive halos (see Eq. Q) and 
to the weak dependence of the peak time, Afpeak, on the halo mass (cf. Fig.|2li. 

The highest luminosity portion of the LFs corresponds to the high tail of the distribution of M, around the mean logM, -logMyir 
relation, at large halo masses. The fit to the data is obtained with a Gaussian distribution with dispersion AlogM, sa 0.3. As 
can be seen in Fig.|3(right panel), for AlogM, « the cutoff in galaxy halo mass (cf. § 2.1) would yield a drastic drop around 
Lb ~ IO'-^'^Lq. We have checked that, for a given mass and virialization redshift of the host halo, a reasonable variation of the 
physical parameters from structure to structure can account for such a scatter. For example, at z ~ 6 and Myir « 2 x IQ^^Mq it 
can be achieved if the clumping factor C varies by a factor of 2, or the seed BH mass varies by a factor of 10, or the strength of 
the QSO feedback eQso varies by a factor of a few. Mahmood et al. (2005) also noted this problem, and empirically solved it by 
inserting a Lorentzian tail in the BH formation rate for halo masses above IO'^'^Mq. 

The observations by Richards et al. (2006, see their Fig. 20) show that the number density of very luminous QSOs with 
^1450 ^ ~27 peaks between z = 2 and z = 3. Our model reproduces this trend. Specifically, the rise from high redshift to z ~ 2.5 
is due to the strong increase of the formation rate of very massive halos, which overwhelms the effect of the decrease of the BH 
mass associated to a given Mvii (see Eq. 13). But at redshifts z < 2 the latter effect dominates, causing the fall in the number 
density of bright sources. 

To reproduce the optical data, we adopt a visibility time Afyis ~ 5 x 10^-10'^ yr, the longer value being favored at redshifts 
z ^ 2. The X-ray LFs require visibility times around Afvis « 3 x 10^ yr at 1-5 < z < 3. The visibility time turns out to be a 
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Fig. 4. — Left panel: predicted QSO LF in tlie optical B-band at redshift 3 for A = 1.7 and different values of the visibility time: Afyis ~ 5 X 10^ yr (best fit, 
solid line), Afyis ~ 10** yr (dashed) and A?vis ~ 3 X 10** yr (dot-dashed). Right panel: effect of the scatter around the mean logM, — logMyir relation. The solid 
curve is the same as in the left panel and con'esponds to a scatter of 0.3 dex, as assumed throughout this paper. The dashed curve is for a scatter of 0.15 dex and 
the dotted one for zero scatter. In both panels data are from Pel (1995, circles), Hunt et al. (2004, triangles), and Richards et al. (2006, squares). 



factor of 5- 10 shorter than the time Afpeak spent by the BH to grow to its final mass. This implies that absorption has a key 
role in determining the duration of the visible QSO phase; in fact, large amounts of dust are present in the pre-QSO phase, as 
witnessed by sub-mm observations. Note that such a phase can also be probed through hard X-ray observations, than can detect 
the emission from the growing supermassive BH. Actually such emission has already been revealed from the center of sub-mm 
selected galaxies (Alexander et al. 2005; Borys et al. 2005); a specific discussion of this issue is presented by Granato et al. 
(2006). 

As stated in § 3.1, the accretion rate onto the central BH has been taken to be at most slightly super-Eddington, X=M,/M^'^'^ < 
4, as expected when it is limited by radiation pressure (see Small & Blandford 1992; King & Pounds 2003). Observational 
evaluations of the Eddington ratio A for high luminosity QSOs with L > 10^^ erg s"' at redshifts z > 2 have been obtained by 
Warner, Hamann, & Dietrich (2004) and Vestergaard (2004). These authors find a significant fraction of QSOs radiating at 
mildly super-Eddington powers. On the other hand, at lower redshifts and luminosities the emission seems to be limited to A < 1 
(McLure & Dunlop 2004; Kollmeier et al. 2006). Note that these results are to be taken cautiously. In fact, they are based on 
the observed FWHM of the Hp, Mgll and CIV broad lines and on an empirical relation linking the size of the broad-line region 
to the QSO luminosity. On top of that, a normalization to the M, - a relation is needed to set the scaling factor in the equation 
yielding the BH mass. This coefficient depends crucially on the kinematic and geometrical properties of the broad line region 
itself (see Onken et al. 2004). Thus the inferred Eddington ratios suffer from large observational and systematic uncertainties, 
which amount to a factor of up to 10 on individual objects, and up to 4 on statistical samples (see Vestergaard & Peterson 2006). 
In our model a variation of the Eddington ratio A with the redshift is necessary to fit the QSOs LFs (see § 6 for a discussion of 
this point). 

If radiation pressure keeps the BH growth at around the Eddington limit, outflows with mass rates ^ Mem are expected, 
as shown in Appendix A. In massive galaxies the outflows remove most of the cold and infalling gas. The model predicts the 
existence of clouds of chemically enriched gas flowing out from the host galaxies (the expelled cold gas). Such outflows have 
indeed been detected by spectroscopic studies of narrow absorption lines associated to QSOs (Srianand & Petitjean 2000; D' 
Odorico et al. 2004). The average metallicity of the gas expelled from large galaxies (Mvir > IO'^Mq) is Z ~ 1 -2Z0, a value 
that after proper dilution contributes to the metal abundance Z ^ Zq /3 of the intergalactic medium in the central regions of 
clusters (see the review by Voit 2004). Note, however, that gas clouds coming from the innermost regions can have strongly 
supersolar metallicities. 

The accreted mass function predicted by our model is plotted in Fig.Q and is found to be in good agreement with the empirical 
determinations by Shankar et al. (2004) and Marconi et al. (2004). This shows that about 60% of the relic supermassive BH mass 
density is in place already at z > 1 .5. The rest of the mass is accreted at later times onto BHs with final mass M, < 10^ Mq. This 
result agrees with observations, showing that at 0.5 < z < 1.5 QSOs of relatively low luminosity, Lx ^3 x 10*'^ erg s"', provide 
about 60% of the X-ray background (Ueda et al. 2003). The 40% in mass accreted at z < 1.5 can account for this fraction if the 
ratio of X-ray to bolometric luminosity increases with decreasing redshift and/or luminosity, as suggested by observations (e.g., 
Wilkes 1994; Vignali et al. 2003; Strateva et al. 2005). The accreted mass function is consistent with the Soltan (1982) argument. 
We have checked this by redoing the calculations of Shankar et al. (2004) with the Eddington ratios adopted in this paper, and 
exploiting the Ueda et al. (2003) LFs corrected for very obscured sources. 

4.2. Host galaxy properties 

The model, interfaced to the GRASIL code, provides also good fits to a variety of galactic properties. Most of the results have 
been presented in previous papers (Granato et al. 2004, Cirasuolo et al. 2005, Silva et al. 2005) and are summarized in Table A2. 
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Fig. 5. — Left panel: The QSO LF in the optical B-band at redshift 2 where A = 1, with the visibility time set to A?vis ~ 10 yr; data points are from Croom et 
al. (2004; circles), from Richards et al. (2005; triangles), and from Richards et al. (2006; squares). Right panel: The QSO LF in the hard X-ray band (2-10 keVs) 
at redshift 2 where A = 1, with the visibility time set to Aiyis ~ 3 X 10* yr; data are from Ueda et al. (2003, shaded area), from Barger et al. (2005, diamonds), 
and from La Franca et al. (2005, arrow). 



Note that the current set of parameters differs from the one adopted in Granato et al. (2004), but the general scenario and in 
particular the results concerning the galaxy properties are essentially unaffected. 

Figure[2shows that small and large halos exhibit quite different behaviors. In halos with masses Myir < IO^^Mq the interplay 
between SFR, BH accretion, and respective feedback processes, limits the growth of the central BH, which never reaches the 
power needed to remove all the gas. Thus the SFR is prolonged over times Af* > Afpeak, i-C-, a substantial amount of stars is 
formed after the peak of the QSO activity. Contrariwise, in large galaxies with Mvh > IO'^Mq the growth of the BH produces 
winds, which are able to stop star formation in a short time w Afpeak, see Fig.|2](right panel). 

On the other hand, in the pre-QSO phase the SFR varies from several hundred to several thousand solar masses per year, 
within halos of mass ranging from 10'^ to lO'^M©. Almost all the associated power is emitted in the far-IR. The model fits the 
850/im counts; in particular, it predicts a large surface density A^(> lmJy,z > 5) w 600 (sq. deg.)"' of pre-QSO host galaxies at 
substantial redshift among bright objects selected at 850 fim (Silva et al. 2005). This large density (compared to that of the QSOs 
at the same redshift) is mainly ascribed to the fact that the pre-QSO phase lasts a factor of 10 longer than the QSO phase, i.e., 
Afpeak/Afvis « 10 at high redshifts. 

A direct consequence of the extremely high SFR in massive halos is the early appearance of galaxies with large mass in stars 
(M* cx SFR X Afpeak at high redshifts). These galaxies, at redshifts corresponding to cosmological times tn > fyh + Afpeak, have 
already gone through the QSO phase and are evolving passively, i.e., without major additions of newly formed stars or large mass 
accretions onto the central BH. 

Note that the 850 /im counts can be reproduced also adopting SFRs much lower than predicted by our model, provided that a 
flatter initial mass function is assumed (Baugh et al. 2005). The latter approach yields a much lower number density of massive, 
high-redshift galaxies. For example, the predicted number density of galaxies with > 10" becomes negligibly small at 
redshift z > 2.5, while in our model it decreases by only a factor of 2 compared to that at z « 1 .5, and it is still significant at z « 5. 
Stronger constraints on this number density would therefore help in discriminating between the two scenarios (e.g., Tecza et al. 
2004; Fontana et al. 2004; Greve et al. 2005; Caputi et al. 2006; Bouwens & Illingworth 2006). 

Our model also naturally accounts for the observed trend of more massive spheroids being richer in metals (see Thomas et al. 
2005). In our model the behavior of the average stellar metallicity (Z^,) as a function of the host halo mass is well described by 
the approximate expression (Z,) « -0.247 -I- 0.023 log (Mv^/Mq). This dependence can be understood as follows. In massive 
galaxies the SN feedback can remove only a small fraction of the cold gas involved in star formation, which is rapidly metal- 
enriched due to the high SFR. When the QSO reaches the peak of its activity, the energy released removes from the halo most of 
both the enriched cold gas and of the infalling metal-poor gas (see also Cox et al. 2006). Afterwards, the gas infall and the star 
formation cease, and the metallicity of the gas does not change anymore. In small galaxies SN explosions are more efficient at 
removing the enriched cold gas but not the infalling metal-poor gas. Thus the former can continuously be replaced by the latter, 
and the metallicity of the cold gas remains low. As a result large galaxies exhibit stellar populations with larger metal content 
than smaller galaxies do. We also stress that the short Af^, « 0.2-0.5 Gyr for Zvh > 4 is responsible for the observed enhancement 
of a elements in massive galaxies (see also Romano et al. 2002; Granato et al. 2004), since the Fe enrichment is strictly related 
to the explosions of Type-la SNae and these are delayed by about 1 Gyr 

The power of the QSO outflows ensures that only about 30% of the initial cosmic baryons are turned into stars in large halos 
(^vir }^ IO'^Mq). For smaller halos the long lasting action of SNae is able to remove an even larger fraction of gas, letting only 
less than 10% of the initial gas to turn into stars (Granato et al. 2004). 

Cirasuolo et al. (2005) found that the observed velocity dispersion function of spheroidal galaxies (Sheth et al. 2003) and 
the Faber-Jackson relation can be reproduced under the hypotheses that the old stellar populations are located in galaxy halos 
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Fig. 6. — Left panel: QSO LF in the optical B-band at z = 1.5, for A = 0.8 and A?vis ^ 10** yr. The data are from Croom et al. (2004; circles), Richards et al. 



(2005; triangles), and Richards et al. (2006; squares). Right panel: QSO LF in the hard X-ray band (2-10 keV) at z ■ 
data are from Ueda et al. (2003, shaded area), Barger et al. (2005, diamonds), and La Franca et al. (2005, arrow). 
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virialized at redshift z > 1.5 and that the relation a « 0.55 Vvh holds; here the velocity dispersion cr refers to the old stellar 
populations in the central galaxy regions. They also showed (cf. their Fig. 4) that the same (T-Vvir relation combined with the 
M, -Mvir relation for redshift z > 1 .5, fits the data in the M, - a plane, with their dispersion (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000; Tremaine et al. 2002). 

Our model predicts a steepening of the relationship in this plane at low masses, since for small halos SNae damp down the SFR, 
thus hampering the growth of the reservoir around the seed BHs and eventually of the BHs themselves (cf. Eq. IA4I : see also 
Granato et al. 2004). This steepening is fully consistent with the upper limits on the BH mass for M33 (Gebhardt et al. 2001) and 
NGC205 (Valluri et al. 2005) obtained by direct dynamical measurements, although the indirect determinations by Greene & Ho 
(2006) suggest a constant slope and a wide scatter at small BH masses. On the other hand, we should also keep in mind that the 
low BH mass portion of the diagram might just reflect the distribution of BH seeds possibly created by merging of smaller BHs 
during the fast accretion epoch and be only weakly affected by the mass accretion. In addition, low BH masses may be easily 
increased by substantial factors during later accretion phases, when host galaxies are disk-(rather than bulge-)dominated. 

5. COMPARISON WITH PREVIOUS MODELS 

Accounting for the wealth of data on the QSO and host galaxy LFs and on their statistical properties and relationships is a 
severe challenge for analytical, semi-analytical and numerical galaxy formation and evolution models. 

The standard methodology adopted by analytical models to estimate the QSO LFs is to convert the DM halo formation rates 
into BH formation rates through a relationship between the mass of the supermassive BH and that of the host halo. Then a 
description of the QSO light curve (usually assuming Eddington-limited accretion) is adopted to derive the QSO luminosity. The 
QSO LFs are then built up according to the approximate expression 



$(L,z) « Afv 



d^A^ST 

dfvir dMvi 



dMv 



dM. 



dM. 



(9) 



The above equation corresponds to Eq. (|5}, under the assumptions that the duty-cycle of QSO activity Afvis is much shorter than 
the cosmological time and that the QSO appears as soon as the host halo virializes, i.e. Afpeak ~ 0. 

Wyithe & Loeb (2003) and Mahmood et al. (2004) estimated the QSO LFs as a function of cosmic time by assuming: (i) no 
delay between virialization and QSO peak luminosity; (ii) a visibility time fvi,s ~ fdi,sc where the dynamical time of the galactic 
disc is fciisc = 0.035 /?vii7Kii ~ 7 x 10^(1 +zT^^^ yr; (iii) halo formation rates derived from the extended Press & Schechter theory 
(Lacey & Cole 1993); (iv) the M, -Myii relation deduced from the relations M, - cr and a-M^^y observed by Ferrarese (2002). 

The first and second assumptions are quite critical at high redshift. If the time scale of QSO activity is associated to the dynam- 
ical timescale of a galactic disc, one has first to wait for the disc to be built up, which implies a delay time At > max[fcooi,fdyn], 
where both times are much larger than fdisc- 

The specification of the DM halo formation rates is also a crucial step for analytical models. Most of the papers in the literature 
adopt the extended Press & Schechter theory (Bond et al. 1991). In this paper we have used instead the rates provided by the 
positive part of the time derivative of the Sheth & Tormen mass function. Note that in both cases the rates are given by the product 
of the mass function (Press & Schechter or Sheth & Tormen, respectively) by a factor, weakly dependent on time (Kitayama & 
Suto 1996). As a matter of fact, mass functions derived from numerical simulations are much better approximated by the Sheth 
& Tormen than by the Press & Schechter theory; the latter at high redshifts and at the large halo masses, relevant for QSOs, can 
be smaller by a factor of 10 (e.g., Springel et al. 2005). 

As for the last assumptions recalled above, Wyithe & Loeb (2003) and Mahmood et al. (2004) assume a M, -Mvh relationship 
steeper and more rapidly evolving than ours. The overall result is that, on average, they predict much more massive BHs at high 
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Fig. 7. — Comparison of the predicted mass function of supermassive BHs associated with massive spheroidal galaxies (solid line) with observational estimates 
(Shankar et al. 2004). The dark shaded area refers to the local supermassive BH mass function considering only early-type galaxies; the light shaded area includes 
also bulges of late-types. 



redshift and for large halo masses (see also the discussion in Shankar et al. 2006). For example, at z = 6 and for Mvir = IO'^^Mq 
we have a BH mass M, « 2 x lO^M© while they have M, « 8 x lO^M©. Recall that we obtain the relation of Eq. (|2j as a 
self-consistent output of the equations in Appendix A, that follow the details of the building up of the central BH mass. 

In summary, the comparison of our results with those of Wyithe & Loeb (2003) and Mahmood et al. (2004) is complex, 
because of quite different assumptions. Despite of these, similar QSO LFs are produced. At high redshift their lower fvis S 10^ 
yr and formation rates of DM halos are compensated by: (i) neglecting the time needed to build up the galactic structure and (ii) 
assuming a relationship M, -Mvii that favors larger BH masses. The difference emerges quite clearly at redshifts around 7-8, 
where we predict a dramatic drop off of the LP at high luminosity, because of the delay (see § 4.1), whereas these authors obtain 
number densities larger than ours by a factor of 10. 

Clear differences are also present on the host galaxy side. In our model the time Afpeak is a very important and active phase, 
within which the massive host galaxies are forming stars at rates around 10^- 10^ Mq yr~' (cf. Fig.Q and the central BH is 
growing. This phase is directly explored by sub-mm surveys and by the observations of the X-ray emission from the growing 
QSOs detected by Alexander et al. (2005). 

More recently, several authors tried to follow the evolution of both QSOs and host galaxies using results from numerical 
simulations, in which the effects of SN explosions and winds from the central active nucleus are taken into account (see Hopkins 
et al. 2006 and references therein). These authors took an approach complementary to that presented in this paper Instead of 
solving equations describing the basic physical processes, they simulated many realizations of collisions between two stable, 
isolated disk galaxies endowed with a central supermassive BH. A two-phase interstellar medium was used to describe star 
formation and SN feedback (Springel & Hernquist 2003). The accretion onto the BHs was estimated from the local gas density 
and sound speed, and limited to the Eddington rate. They assumed that 5% of the bolometric luminosity derived from accretion is 
transferred to the surrounding gas, as thermal energy. They found that during a major merger the central BHs are fed with enough 
gas to yield a luminosity depending exponentially on time. Their approach includes also estimates of absorption by dust. The 
energy injected by the QSO in the interstellar medium is enough to unbind the gas itself on a time scale Afvis ~ 10^ yr, within 
which the luminosity varies from about lO""* Lpeak to a maximum Lpeak, with a weak dependence on Lpeak [cf. Eq. (7) of Hopkins 
et al. 2006]. Note that these simulations do recover the transition from obscured inflow/BH growth to QSO-driven blowout, as 
we have simply modeled from the relevant accretion/feedback equations. 

The time scale Afyis for the growth of the BHs derived from numerical simulations of low-redshift galaxy mergers at large halo 
mass has the same meaning as our time A/peak when the maximum luminosity is attained. However, for large masses and high 
redshifts the two timescales have quite different values, AfvL, « 10^ yr > Afpeak ~ 0.2-0.5 Gyr forzvii > 4. The long time before 
the peak luminosity is reached has strong implications at high redshifts. In particular, Afvis ~ tn for z > 4. Thus the model of two 
colliding discs takes too long to form supermassive BHs at high z (recall that the discs themselves need at least a time /?vir/Kir 
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to set up). Therefore the extension of the model at z > 2 is quite complex, as discussed by Hopkins et al. (2005) in their § 3.2. 
Moreover, it is likely that several properties inferred from this kind of low-redshift merger simulations cannot be extrapolated to 
high redshifts, where the very idea of binary galaxy mergers likely breaks down. However, it is interesting to note that at lower 
redshifts, z < 1.5, the two times Afvis and Afpeak are comparable. Our model is not aimed at following the evolution of QSOs at 
z < 1. This later phase is better described by models that include 'interactions' (minor mergers, fly-bys, disk instabiUties, etc.) 
among galaxies as triggers of nuclear activity (e.g., Kauffmann & Haehnelt 2000; Cavaliere & Vittorini 2000; Menci et al. 2003). 

6. DISCUSSION AND CONCLUSIONS 

In their thoughtful review paper, Brandt & Hasinger (2005) list some 'key outstanding problems' of AGN astrophysics that 
need investigation. These include: the detailed cosmic history of supermassive BH accretion; the nature of AGN activity in 
young, forming galaxies, and the connection between supermassive BH growth and star formation in sub-mm galaxies. 

In this paper we address these issues in the framework of the anti-hierarchical baryon collapse (ABC) scenario developed by 
Granato et al. (2004). We have shown that the condensation of baryons within DM halos both in stars and in BHs can be described 
with a simple physical model, whose main equations are listed in Appendix A. The model yields the time dependence of the SFR, 
of the accretion rate M, onto the central BH, of the total mass in stars and of the final BH mass M,, for any given halo mass 
and virialization epoch. The GRASIL code then provides the SED from X-ray to radio bands of evolving stellar populations 
(successfully reproducing the epoch-dependent galaxy LFs in different spectral bands, as well as a variety of relationships among 
photometric, dynamical and chemical properties, as shown in previous papers; see Table A2), while observationally determined 
bolometric corrections allow us to convert the accretion rates onto supermassive BHs into luminosities in optical and X-ray bands. 

In the ABC scenario, the growth rate of BHs is proportional to the SFR, and the latter is more effectively slowed down by 
SN feedback in smaller halos, resulting in a relatively more efficient growth of more massive BHs, especially at higher redshifts. 
Both star formation and BH growth in massive halos are stopped by the feedback from the active nucleus as soon as it becomes 
powerful enough to sweep out the residual interstellar medium. 

The BH accretion occurs in host galaxies with very intense star formation, and is heavily dust obscured for the majority of 
its duration; however, most of the final BH mass is gained during the last 1-2 e-folding times in unobscured conditions, as 
demanded e.g. by the Soltan (1982) argument and by the analysis of Yu & Tremaine (2002). For massive objects, the growth 
phase lasts 15-20 e-folding times of Eddington-limited accretion, but the AGNs are detectable by current hard X-ray surveys 
only in the last several e-folding times, i.e. have hard X-ray visibility times of ~ 3 x 10*^ yr. Interestingly, Borys et al. (2005) 
estimate that the BH masses associated to the X-ray emitting AGNs detected in sub-mm galaxies are, on average, ~ 50 times 
lower than those associated to local spheroidal galaxies with similar stellar masses, as expected if they are about 4 e-folding times 
before the maximum, for Eddington-limited accretion. At still earlier times, sub-mm galaxies are expected to host intrinsically 
weak and highly obscured nuclei, undetectable with current X-ray telescopes. 
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The model accounts for the hard X-ray AGN LFs at various redshifts. A crucial ingredient, to this end, is the rapid quenching 
of the accretion rate onto the biggest BHs, predicted by the model as the consequence of the AGN feedback which sweeps out 
the interstellar medium. 

The optical (B-band) visibility time is shorter than that for the more penetrating hard X-rays. The observed epoch-dependent 
B-band LFs are accurately reproduced for a Afvis of order of the e-folding time, indicating that only when the AGN is approaching 
its maximum luminosity it can clear up the surrounding region. The redshift dependence of the space density of optically bright 
QSOs is controlled, in the ABC scenario, by two competing factors. On one side, according to the hierarchical clustering 
paradigm, the formation rate of very massive halos hosting them is increasing rapidly with decreasing redshift. On the other 
side, the BH to host halo mass ratio decreases with decreasing redshift, so that at low redshifts relatively more massive hosts 
are required for a given BH mass. But the galactic halo mass function sinks down exponentially for large halo masses, and is 
actually cut off at Mvir « 2 x IO'^'Mq. In our model, the first factor dominates for z > 2.5 and the second dominates at lower z, 
thus accounting for the increase in the bright QSO space density with decreasing z down to z « 2.5, as well as for the subsequent 
decrease. 

The redshift dependence of the bolometric QSO LP is illustrated by Fig.|8l It is broadly reminiscent of luminosity evolution, 
but with significant deviations from that simple description particularly at high luminosities. We note in particular the flattening at 
high luminosity around z w 5, borne out by SDSS data (Fan et al. 2001) and the very sharp drop (by about 2 orders of magnitude) 
of the QSO LF between z « 6 and z « 8, due to the dearth of massive halos at high redshifts. 

As stated in § 3.1 and 4.1 our model requires a mild variation of the Eddington ratio A with redshift. This may be a matter 
of concern, since the value of A is determined mainly by the environmental conditions very close to the supermassive BHs. The 
main motivation for resorting to super-Eddington accretion at high z is to account for the observed space density of very luminous 
QSOs at z > 5 (an acceptable fit to lower redshift data can be obtained keeping Awl). Since, as indicated by Eq. (jSj, QSOs seen 
at the cosmic time t are associated to halos virializing at the earlier time fvh ^t- Afpeak, and the density of such halos drops with 
increasing Afpeak at high z, to have a sufficient number of massive enough BHs at z > 5 we need them to grow faster than implied 
by Eddington-limited accretion. The values of A mentioned in § 3.1 provide the necessary shortening of Afpeak and also give 
high enough luminosities without resorting to too large BH masses. On the other hand, if we want to preserve the relationships 
between the properties of the QSO and of the host galaxy, Afpeak must be long enough to allow the formation of a sufficiently 
massive stellar component. 

The redshift dependence of A may thus be seen as a way to parameterize aspects of the physics of the accretion/emission 
process not properly taken into account by our simple model. For example, a major improvement could be obtained by taking 
into account the results of the photon trapping theory (e.g., Blandford 2004). According to the latter, radiative emissions close to 
Eddington can be originated by strongly super-Eddington accretion rates (e.g., those allowed in our model by viscous accretion). 
This process seems promising because it may naturally allow short delay times at high redshift even with an Eddington-limited 
radiative output. However, details of the mechanism are still uncertain and additional work is needed to quantitatively include it 
in our model. Non-radiative accretion may also speed up the BH growth. On top of that, the values of the Eddington ratio could 
have a substantial scatter due to different environmental conditions around the supermassive BHs. As a result, the samples of 
high redshift QSOs may be biased towards higher Eddington ratios. 

We note that since luminous high-z QSOs are associated to very massive halos, the model implies that their clustering properties 
are similar to those of massive spheroidal galaxies and of the bright sub-mm galaxies detected by the SCUBA surveys, consistent 
with the results by Porciani et al. (2004) and by Croom et al. (2005). This issue will be investigated in a future paper 

During the relatively short phase when massive galaxies are in the process of expelling their gas and dust, their nuclei should 
appear as powerful X-ray sources, Lx > 10"^^ erg s~', but still somewhat obscured in the optical band; they thus may resemble 
Type 2 QSOs. 

Spectroscopic studies, aimed at determining velocity, amount and chemical abundances of the gas ejected during the QSO 
phase, would be extremely informative. In fact, the model predicts that at early cosmic times large amounts of metal-enriched 
gas were evacuated from the galaxy halos to the intergalactic medium. 

We warmly thank S. Bressan, A. Celotti, and M. Cirasuolo for helpful discussions, and the anonymous referee for a careful and 
perceptive report, rich of illuminating comments. This work is supported by grants from ASI, INAF and MIUR. 
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APPENDIX 

BARYONIC PHYSICS IN OUR GALAXY EVOLUTION MODEL 

In this Appendix we provide an overview of the physical model by Granato et al. (2004), recalling the basic equations and 
parameters that control the evolution of the baryonic component. 

The baryonic content of a given DM halo with mass Mvh is partitioned in three gaseous phases: a hot diffuse medium with 
mass Minf infalling and/or cooling toward the center; cold gas with mass Mcoid condensing into stars; low-angular momentum gas 
with mass Mres stored in a reservoir around the central supermassive BH, and eventually viscously accreting onto it. In addition, 
two condensed phases are present, namely, stars with a total mass and the BH, with mass M,. 

The evolution of the baryonic content is described by the system of differential equations 

Minf = -Mcond-Mi^r , 

Mcold=Mcond-M*-Mr.s-MZi-M^M - (Al) 
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TABLE Al 
Model Parameters 



Parameter 


Value 


Short description 


c 


7 


clumping factor 




0.05 


strength of SN feedback 


OIRD 


2.5 


strength of radiation drag 


TO 


1 


zero-point of optical depth 




lO^M© 


mass of BH seed 


^accr 


10-2 


strength of viscous accretion 


A 


0.8-4 


Eddington ratio* 


V 


0.15 


radiative efficiency 


eQso 


1.3 


strength of QSO feedback 



Note. — * We let the maximum allowed Eddington ratio A to depend on the redshift as: A = 4 for z > 6, A = 3 for 5<z<6, A = 1.7for3<i;<5, A=l for 
2 <z < 3, and A = 0.8 for 1.5 <z < 2; the empirical fit A(z) f5i -1.15 + 0.75 (1 +z) works as well in the^redshift range h5<z<6. ~ ~ 



Mres=Minflow-M. , 

where Mcond is defined by Eq. (A2). At the virialization redshift Zyir we set Mjnf « Mvir/6 corresponding to the cosmic baryon to 
DM ratio and featuring a primordial chemical composition, M, « Ml''''^, M'^^'^'^ being a BH seed mass originated by some process 
in the early universe, and the other baryonic components to zero. 
The hot gas condenses/cools down at the rate 

max[fcool,fdyn] 

determined by the dynamical and cooling timescales fdyn and fcooi at the virial radius. The latter includes the appropriate cooling 
function (Sutherland & Dopita 1993) and allows for a clumping factor C in the baryonic component. Since at high redshifts 
major mergers between massive halos are very frequent, we neglect here the effect of the angular momentum, since it is lost by 
dynamical friction through mergers of mass clouds Mc on time scale t^p « 0.2(^/ln^)fdyn, where f = M^-„/Mc (see e.g. Mo & 
Mao 2004); major mergers imply ^ ~ a few. In this context it is also worth noticing that the gas that collapses and cools enough 
to form stars is only a fraction of less than 30% of that associated to the virialized halo (cf. Fig.[3. 
The cold gas turns into stars at the rate 

/ u 1 ' (A3) 

J max[fcool,fdyn] 

where now fcooi and fdyn refer to cooling and dynamical time, respectively, in the radial mass shell of mass dMcoid- 
The energy feedback from Type-II SNae removes gas from the cold phase at the rate 

where the efficiency esn is a parameter of the model. The amount of cold mass removed is proportional to the number of SN 
explosions (assuming an average energy of around 10^' ergs per SN) hence to the SFR, and inversely proportional to the depth 
of the halo potential well. 

The stellar radiation drag (see Kawakatu & Umemura 2002) triggers the inflow of cold gas into a reservoir of low-angular 
momentum around the central supermassive BH, at the rate 

Minflow ~ auD X 10"^M*(1 -e"^) M© yr"' , (A5) 



where the effective optical depth is given by 

The metallicity Z(f ) is computed from the code self-consistently, while the coefficients a/to and tq are parameters (see Granato 
et al. (2004) for further details). 

The viscous time fyisc = Rccrit fdyn^'^'^'* ^'^^ '■^^ accretion from the reservoir to the BH is related to the dynamical time fdyn^"'^ °f 
system BH -I- reservoir, and to the critical Reynolds number Rccrit ^ 10^- 10^ (see Granato et al. 2004 for details). Thus the gas 
in the reservoir accretes onto the supermassive BH at the rate 

Mr «5x 10^ fc_f^^Vf^V'7i+—V'' ^oy-' ■ (A7) 

VSOOkms-V \M,I \ M,,J °^ 
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TABLE A2 
Overview of Model Results 



Property or Statistics 



External Inputs* 



Reference 



QSOsandBHs 

M, -Mvir 



A?, 



peak 



Opt. LFs (z) 
X-ray LFs (z) 
BH mass function 
M.-o- 



no add 
no add 
A?vis fa 

A?vis Ri 

no add 
no add 



5 X 10'' yr, AlogM, ; 
3 X 10* yr 



Fig. 2 (left) 
Fig. 2 (right) 
Figs. 3-6 
Figs. 5-6 (right) 
Fig. 7 

Fig. 6of G04 



Spheroidal Galaxies 
VDF 

Faber-Jackson 
LF(z) 

Metal abundances 
K-band counts 
850 /um counts 
EROs 



<T/Vvir 

IMF, SEDs (GRASIL) 
IMF, chemical yields 
IMF, dust modeling (GRASn.) 
IMF, dust modeling (GRASn.) 
IMF, dust modeling (GRASn.) 



Fig. 1 of C05 
Fig. 2ofC05 
Fig. 10ofG04 
Fig. 9ofG04 
Fig. lofS05 
Fig. 12ofG04 
Fig.6ofS05 



Note. — * Input quantities not included in Table Al. G04; Granato et al. (2004), C05: Cirasuolo et al. (2005), S05: Silva et al. (2005). 



Again, the coefficient A:accr is a model parameter. For the plausible range of 10~^ < ^accr ^ 10"^, the accretion rate Mi'*"^ easily 
exceeds the Eddington rate for high redshifts and massive halos. In this situation we may have: (i) very rapid growth of the 
central BH in a dynamical timescale, followed by a slower accretion when less mass is available in the reservoir; (ii) BH growth 
limited by radiation pressure, followed by strong mass outflows M„ ~ (e.g., Small & Blandford 1992; King & Pounds 2003; 
Begelmaim 2004). We select the latter option and the actual accretion rate is set to 



M. =min[M™^AM^''''] , 



where 



1.3 



lO^Mf: 



Mq yr 



(A8) 



(A9) 



is the Eddington rate and A is the maximum allowed Eddington ratio. We let the maximum allowed Eddington ratio A to depend 
ontheredshiftas: A = 4 for z > 6, A = 3 for 5 < z < 6, A = 1.7 for 3 < z < 5, A = 1 for 2 < z < 3, and A = 0.8 for 1.5 < z < 2. The 
empirical fit A(z) w -1.15 + 0.75(1 +z) works as well in the redshift range 1-5 < z < 6, see also Table Al. 

In the presence of Eddington-limited or super-Eddington accretion, QSO-driven outflows are expected in the form of winds 
that can form just above the disc by a combination of radiation and gas pressure. Large mass outflows have been confirmed 
by X-ray observations of BAL QSOs (Brandt & Gallagher 2000; Chartas et al. 2003). Following the model by Murray et 
al. (1995) the asymptotic speed and the mass outflow rates are estimated as v,v w O.O6c(r„/lO"'cm)-'/2(M,/M0yr-i)i/2 and 
M„ « 2.4 fcN22 (r,v/lO^^cm)'''^(M,/M0yr"^)'/^ Mq yr"\ respectively; here c is the speed of light, fc is the covering factor, 
N22 is the hydrogen column density normalized to 10^^ cm~^, and r„ is the launching radius. Then the kinetic luminosity of the 
outflow reads 



i 2.5 X 10 



,44 



(fc_\ (N22\ ( r„ y'/^ / M. A 
VO.l y V 10 / VlOi6cm7 \MQ/yrJ 



3/2 



Eventually, QSO outflows remove gas from both the hot and the cold phases at the rates 



ergs 



3/2 



Mm. 



.cold 



Mg/yry Mcoid+Minf 



Mq yr 



(AlO) 



(All) 



with eQso = (//i/0.5)(/c/0.1)(A^22/10), where //, is the fraction of the kinetic luminosity that couples to the gas. Note that most of 
the power L/j = fyLg is transferred to the gas near the peak of the BH accretion rate curve. Such a power is only a small fraction 
of the Eddington luminosity (L^j/LEdd) « 1-4 x lO"^e2so(M,/lO^M0). 

The basic outputs of the model are the evolution of the SFR M*(f) and the BH accretion rate M,(t) as a function of the galactic 
age t for any given value of the halo mass Mvii and of the virialization redshift Zvir- Once the star formation history of a galaxy 
has been computed, its luminosity in any chosen band is obtained as a function of / from the GRASIL code, which yields the 
chemical and the spectro-photometric evolution from the radio to the X-ray band, allowing for the effect of dust absorption and 
reradiation. 



